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Abstract. 

We demonstrate that a phonon stopband can be synthesized from an aperiodic 
structure comprising a discrete set of phononic filter stages. Each element of the set has 
a dispersion relation that defines a complete bandgap when calculated under a Bloch 
boundary condition. Hence, the effective stopband width in an aperiodic phononic Hlter 
(PnF) may readily exceed that of a phononic crystal with a single lattice constant or a 
coherence scale. With simulations of multi-moded phononic waveguides, we discuss the 
effects of finite geometry and mode-converting junctions on the phonon transmission 
in PnFs. The principles described may be utilized to form a wide stopband in acoustic 
and surface wave media. Relative to the quantum of thermal conductance for a uniform 
mesoscopic beam, a PnF with a stopband covering 1.6-10.4 GHz is estimated to reduce 
the thermal conductance by an order of magnitude at 75mK. 
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The ability to taylor the phonon dispersion curve in a meta-material phononic crystal 
(PnC) is ideal for managing heat transport. In addition to the phonon group velocity 
that can be engineered in a PnC to control the phonon flux [T], the existence of 
complete phonon bandgaps promises a new class of very low thermal conductance meta¬ 
materials ISli- In general, topology optimization is employed to maximize the bandgap 
width, Au, arising from a unit cell arranged in an inhnite periodic lattice. This approach 
is inherently limited by the fact that the lattice constant, a, hxes the length scale for 
resonant scattering and Bragg interference in the PnC. In addition, as the gap center 
frequency, ujq, is increased, the algorithms tend to generate progressively more complex 
shapes with feature size less than the lattice constant a [H El E]- Guided by these 
observations, we explore the hltering properties of a set of discrete coherence length 
scales combined to from a phononic hlter (PnF) as illustrated in £gure[^ In the periodic 
limit, the phonon dispersion of the geometry in each hlter stage exhibits one or more 
complete phonon bandgaps. The bandgaps can be thought of as a hlter pole, arising 
either from Bragg interference or resonance within the phononic structure. Thus, the 
ehective stopband of a PnF is a consequence of several hlter poles distributed across a 
wide bandwidth. 

In electromagnetism, the use of hnite-sized structures as hltering elements is a 
concept that is well studied, for example, in the form of waveguide chokes, stepped 



Figure 1. Example one-dimensional aperiodic phononic filters with wide stopbands. 
The individual filter stages employed are selected from geometries that result in 
complete phonon bandgaps when simulated in the periodic limit. The total stopband 
of the SxN stage filter extends from the lowest bandgap edge of cell A, to the 
highest bandgap edge of cell D. Each filter realization is utilized to explore a particular 
characteristic of a PnF, including a wide stopband in (a), and effect of mode conversion 
on low frequency phonons in (b). The inset shows a perspective view of a stepped-width 
unit cell. The thickness t is 300nm and constant across the planar phononic structures. 
The lengths A, B, B’, C, and D, are 780, 645, 570, 550, and 400 nm respectively. The 
width a is equal to the length in each stage. The uniform input/output beams in each 
filter are 300x300 nm^ in cross section, and matched to perfectly absorbing layers in 
finite-element computations. The dashed arrow indicates the transmission sense -\-k of 
the mechanical modes as defined by waveports in the model. 
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impedance filters [7j, and photonic structures [8]. A simple example is provided by 
a lumped capacitor and inductor combined to form a notch hlter stopband. Adding 
additional lumped or distributed circuit elements to the filter can be used to adjust 
the response to achieve the desired stop-band specification. The work described here 
can be viewed in this context as the acoustic analog of a multi-mode transmission line 
structure. We refer to Brillouin ^ for a detailed discussion of how elastic bodies can be 
thought of or decomposed into elastic reactive elements. In contrast to electromagnetic 
single-mode waveguides, phonon waveguides are inherently multi-moded with at least 
four acoustic modes for a ID beam (out-of-plane flexural, in-plane flexural, torsion, and 
compression) [TU], where D is the dimensionality of the reciprocal space. In practice, 
both the acoustic and higher order (optical) modes need to be considered. 

This work is motivated by the need for ultra-sensitive and fast cryogenic thermal 
sensors for future space-based far-infrared telescopes. The thermal-fluctuation-hmited 
noise-equivalent-power of a bolometer is given by AksT^G^ where Tq is the operating 
temperature of the bolometer, and G is the thermal conductance of the bolometer to 
the cryogenic environment m- For state-of-the-art refrigeration systems capable of 
cooling a kilo-pixel detector array to 100 mK, reducing G provides a path for reaching 
~ 10“^® W/v^Hz sensitivity without increasing cryogenic complexity. A wide stopband 
PnF can potentially limit the heat transport to the desired level in a compact structure 
less than 10 pm in lateral dimensions, an advantage given that a focal plane can be 
packed efficiently, and the hlters will contribute negligible thermal mass to the sensor. 

In Section we describe the framework of multi-moded ID low-pass phononic 
filters. We solve the isotropic elastic wave equations using the finite-element method 
(see Appendix). In Section]^ we assess the effect of finite-geometry, and compute the 
transmission spectrum of the four lowest acoustic modes of a uniform beam through a 
phononic filter as illustrated in figure [^a). In Section]^ we describe the effects of a 
mode-converting junction (figure[^b)) on the transmission of the compliant out-of-plane 
flexural and compression modes. In Section we discuss application of the method in 
terms of thermal conductance and sensitivity of a bolometer thermally isolated with ID 
PnFs. 

2. PnF Stopband 

In a periodic phononic crystal, the bandgap center frequency, ojq, scales inversely with 
the size of the unit cell. This fundamental relationship that arises from the coherent 
properties of the crystal is the basis of the phononic hlter approach. However, for 
a planar ID structure with constant thickness, merely scaling the lattice constant, a, 
does not guarantee a precise shift in ojq. The frequency response is a function of the 
lattice constant to thickness ratio [3]. Thus, we simulate phononic crystals with a 
Bloch boundary condition, select geometries that exhibit complete phonon bandgaps, 
and repeat the calculations as a function of a. Figure shows that for the material 
and unit cell geometry explored here, bandgaps across a continuous bandwidth from 
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1.6-10.4 GHz may be identified. Although other shapes could be used, the choice of a 
planar stepped-width unit cell reflects the relatively large complete phonon bandgaps 
that can be produced in a periodic crystal of this from. The step in width corresponds 
to a change in mechanical driving point impedance [iniiig. 

The unit cell geometries selected in figure can be cascaded to form a finite-sized 
phononic filter, as shown in figure [T] The filter would exhibit a phonon stopband that 
is synthesized from numerous poles distributed across a wide bandwidth. Poles arising 
from resonant scattering or Bragg interference may contribute to the response. We 
note that the stopband would formally transition to a complete phonon bandgap if a 
phononic filter is repeated as infinite periodic array. 

For the management of heat transport, the number of filter stages in a PnF can be 
treated as a free parameter when thermal isolation is the only objective, as is the case in 
thermo-electric generation. In applications where energy storage is also important, the 
rejection level and width of the stopband must be traded against increased thermal mass 
and parasitic effects, such as enhanced acoustic loss due to coupling to evanescent modes 
and fabrication imperfections. Thus, in the approach described here, SxN is minimized 
and the overlap of the bandgaps is optimized to reduce the thermal occupation of 
standing-wave modes where doj/dk ~ 0, and inter-band leakage. 

The fractional bandwidth of the synthesized stopband in an aperiodic PnF can 
readily exceed 100%. To evaluate the usefulness of a large rejection bandwidth, consider 
the broad phonon emission spectrum of a blackbody source through a ID phononic 
channel, as shown in figure For maximum effectiveness, the filter stopband is biased 
towards the Rayleigh-Jeans limit. However, even at 100 mK, 10% of the thermal power 
is carried by phonons with frequencies greater than 10 GHz. Hence the typical 30% [ 13 ] 
fractional bandwidth in a phononic crystal will have limited impact on the thermal 
conductance in a mesoscopic structure [3|. 

Above ID, the PnF stopband should be placed over the peak frequency ujth of the 
Planck distribution. In particular, the normalized bandwidth Au/uth must be greater 
than unity to achieve substantial reduction in thermal conductance. To achieve this goal 
in a coherent structure such as a PnF, the cell size should span the appropriate range to 
enable a wide stopband. In a uniform beam with boundary-limited scattering [TT], the 
spectral density of the surface features needs to be sufficiently large in the waveband 
of interest to ensure diffusive phonon propagation. Fabrication tolerance ultimately 
determines the physical limit (coherent or diffuse) that can be successfully accessed to 
control heat flow. 

In a coherent phononic filter fabricated using electron-beam lithography, the degree 
of surface asperity is expected to induce negligible phonon decoherence. It has been 
shown that in electromagnetic meta-materials, decoherence becomes important when 
the scale of the surface roughness approaches the wavelength of the propagating 
electromagnetic mode [T5|. A similar argument readily applies to elastic meta-material 
waveguides. Table summarizes the relevant roughness parameters in nano-meter 
features fabricated with electron-beam lithography [16]. The root-mean-square and 
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Figure 2. Complete bandgaps of a stepped phononic cell geometry calculated 
under Bloch boundary condition (lower panel). A cascaded aperiodic or quasi-periodic 
structure formed from a combination of these cells will exhibit a stopband covering 
1.6 to 10.4 GHz (grey region). The cell dimensions are within the fabrication limits of 
state-of-the-art electron-beam lithography. A, B, C, and D refer to the cell geometries 
of the PnF shown in figure [^a). The beam width w is 50 nm in all cases. The 
stopband is compared to the effective bandwidth of a one-dimensional thermal source 
(upper panel). g{i',T)/gQ represents the fractional cumulative thermal conductance of 
an acoustic phonon mode, where g{v,T) = J^[dP{v^T)/dT\dv^ P{v,T) is the one¬ 
dimensional phonon spectral density, and go = g{0,T). 


correlation length in beam edge roughness are subdominant to the thermal wavelength 
at 100 mK by two orders of magnitude, Xth ~ 27ihv/kBT, where v is the typical speed 
of sound of out-of-plane flexural modes in a dielectric (~5000 ms“^). Hence, cooling 
phononic devices to sub-Kelvin temperature offers the opportunity to fully test and 
measure the coherent long-wavelength meta-material properties of these structures at 
the mesoscopic limit [171 [3l HE] • 

3. Finite Geometry Effects 

The Bloch boundary condition is convenient for computing the elastic eigen-modes of 
an inhnite crystal, however, a phononic crystal is truncated in a practical realization. 
The effect of hnite geometry, or number of hlter stages on the phonon transmission 
coefficient, is of therefore particular interest. 

Figure [^illustrates the typical change in transmission of the compliant out-of-plane 


Table 1. Typical surface roughness characteristics induced by electron-beam 
lithography m- For reference, the phonon thermal wavelength at 100 mK is shown. 


Surface Roughness RMS 

(Ts ~ 5 nm 

Edge Roughness RMS 

(Je ~ 5 nm 

Edge Roughness Correlation Length 

~ 100 nm 

Thermal Wavelength 

Xth ~ 2000 nm 
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flexural and compression modes of a uniform beam through a quasi-periodic A^-stage 
hlter (5 = 1). As N is increased, the filter response asymptotically approaches that of 
an infinite array: the upper and lower stopband edges shift towards the ideal bandgap 
edges, and the transmission nulls deepen. Thus, as the crystal is truncated to a finite 
size, the complete phononic bandgap transitions to a stopband in which the phonon 
transmission of the propagating acoustic modes are reduced. When the wavelength is a 
signihcant fraction of the hlter stage size, the minimum of the transmission coefficient in 
the bandgap regions scales exponentially with increasing N, as expected of evanescent 
decay in a coherent hlter [ZlIISKla]. 

The required rejection of a phononic hlter is application dependent, and typically 
determined by the signal strength, sensor noise and bandwidth. Three to hve hlter stages 
may be sufficient to provide the level of rejection needed for many applications [13]. As 
described in section in a practical implementation, increasing the number of hlter 
stages should be traded against increasing thermal mass and decoherence arising from 
geometric imperfections. 

To illustrate a phonon stopband within a hnite aperiodic geometry, we examine 
the hlter shape shown in hgurej^a), with 5=4 and N=2. 5 and N ehectively dehne 
the stopband width and the level of rejection respectively. The transmission coefficient 
of the four acoustic modes is shown in hgure The lower stopband edge for the out- 
of-plane hexural and compression modes occurs at 1.6 GHz, the lowest frequency of 
the hrst phonon bandgap in the hlter calculated in the periodic limit (see hgure [^. 
The in-plane hexural and torsion modes are strongly scattered due to the quadratic 
dependence of the mode impedance on the beam width [201 E|. The out-of-plane 
hexural mode impedance also depends quadratically on thickness, however, because of 
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Frequency (GHz) Frequency (GHz) 


Figure 3. Power transmission of the compliant acoustic modes through a quasi- 
periodic filter with N hlter stages (5=1). Referring to the unit cell geometry in the 
inset of hgure [2| a = 780 nm, I = 75 nm, and w = 50 nm. The grey region shows the 
bandgaps of the phononic crystal when calculated in the periodic limit. 
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Figure 4. Power transmission through a PnF with the geometry shown in figurej^a). 
The grey region is formed by the overlapping bandgaps of each hlter stage calculated 
in the periodic limit (see figure . The solid line at -30 dB is for reference. The 
frequency resolution in the simulation is 17 MHz. 


the constant thickness of the geometry, these modes experience less rejection by the 
hlter. The impedance of compression modes is linear in thickness and width. Thus 
out-of-plane hexural and compression modes are most susceptible to parasitic effects in 
the phononic hlter geometry described here, although even at N=2, the transmission 
of the low-pass hlter is below -30 dB from ~2 to 8 GHz. We emphasize that the hnite- 
element computations take into account the rich mode coupling between all the modes 
supported by the PnF, including evanescent modes. 

4. A Mode Converting Junction 

Below the lowest stopband edge frequency of the PnFs, the compliant out-of-plane 
hexural and compression modes have near unity transmission. For the out-of-plane 
hexural mode, the transmission may be ehectively limited if the thickness is varied 
across the hlter geometry |2U]. In applications where the device thickness is uniform, 
however, mode-conversion to the torsional and in-plane hexural modes may be employed 
to reduce the transmission of both compliant modes. 

An example PnF with a mode-converting junction is shown in hgurej^b), where 
the 90° bend scatters out-of-plane hexural into torsional modes, and compression into 
in-plane hexural modes. By symmetry, the reverse scattering processes can also occur, 
however, the low transmission of the in-plane hexural and torsional modes across the 
PnF transfers negligible power. Relative to a PnF with the same hlter stage dimensions 
but a straight prohle, the 90° bend reduces the sub-Gigahertz average transmission by 
a factor of 5.5 for the out-of-plane hexural mode, and a factor of 4 for the compression 
mode, as shown in hgure 
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Figure 5. Power transmission through the phononic filters shown in figure a) and 
(b), solid and dashed lines respectively. S = 3 and TV = 2 in both filters. The grey 
region is the bandgap of the first filter stage (a = 780 nm) calculated in the periodic 
limit (see figure]^. 


5. Thermal Conductance of a PnF 


At 100 mK, the thermal wavelength exceeds the cross-sectional dimensions of the 
waveguides in hgure Thus, the thermal conductance of the uniform beam through 
the PnF can be estimated by incoherently summing over the contribution from the 
quantized ID propagating modes, 

exp (a; 




Tr,(x) 


;dx, 


( 1 ) 


(exp(a;) — 1)^ 

where x = hu/ksT, = x’^ksT/U is the cut-off angular frequency of mode a, and Ta is 
the power transmission coefficient through the PnF. For the four lowest acoustic modes, 
a = [1,4], = 0. These modes contribute to most of the thermal conductance when 

T-ChcaJ/fc^. Equation]^ is valid in the absence of acoustic loss, or thermalization, in a 
PnF. This approximation is reasonable given the expected level of decoherence induced 
by geometric imperfections in a practical device, as discussed in section The upper 
bound on the total ballistic quantized thermal conductance, Gq, can be evaluated with 
Ta (a;)=l, Va. At 100 mK, the thermal conductance fraction carried by an acoustic mode 
extends far above the Rayleigh-Jeans equivalent bandwidth, 'xksT/{12%) = 3.4 GHz, as 
illustrated in hgure 

To establish a hgure of merit, we use equation [T] to calculate the thermal 
conductance of a uniform beam with and without a PnF. Also of interest is the hltering 
performance of a PnF relative to a phononic crystal with a hxed lattice constant. For 
the latter, the cell geometry {a,l} = {780,75} /rm in hgure is chosen as an ideal 
representative case with an inhnite number of hlter poles leading to unity transmission 
outside and zero transmission within the complete bandgaps. The fractional widths are 
20% for the hrst, and ~10% for the second and third bandgaps respectively. 
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Figure 6. (a) Thermal conductance of various ID phononic waveguides, normalized 

to go = kgTrT/{6h). The uniform beam cross-section is 300x300 nm^. All waveguides 
are matched to thermal baths at both ends. a;5/27r=8.5 GHz is the cut-off frequency 
of the first waveguide mode of the uniform beam. A = 4 and = 2 for the straight 
and bent phononic filters, (b) The thermal-fiuctuation-limited noise-equivalent-power 
of a bolometer suspended by three filtered beams with conductance given in the left 
plot. An NEP < 4 X 10“^® W/-\/Hz is achievable at lOOmK. 


Limited by the mesh size of the hnite-element computations, we approximate 
the phonon transmission in the maximum possible synthesized stopband of the PnFs 
explored here (1.6-10.4 GHz, see hgure as follows. We compute the phonon 
transmission coefficient through the hlter prohles in hgure with 5 = 4 and N = 2 
up to 8.5 GHz. The phonon transmission for every mode of the uniform beam in the 
8.5-10.4 GHz range is set to -30 dB, a reasonable assumption given the expected rejection 
in each hlter is an order of magnitude below this level. Above 10.4 GHz, the transmission 
is set to unity. The contribution of the optical waveguide modes of the 300x300 nm^ 
beam up to 30 GHz is included in all calculations. 

Relative to the total quantum of conductance Gq for the uniform beam, an order 
of magnitude reduction in thermal conductance is achievable with a PnF below 75 mK, 
as shown in hgure |^a). In absolute terms, the estimated thermal conductance of a PnF 
with a 90° bend is 94 and 27 fW/K at 100 and 75 mK respectively. 

Depending on the application, a practical phononic hlter may have one or both 
ends attached to a much wider membrane, essentially forming a discontinuous junction. 
At the low-frequency (long-wavelength) limit, the transmission of the acoustic modes 
tends to zero as a power law, I2D1. p is 1/2 for the out-of-plane hexural mode, and 
1 for the compression mode. Thus, the sub-Kelvin thermal conductance in hgure [^a) 
may be considered as an upper bound on the performance of a PnF as modeled with 
perfectly matched thermal baths at both ends of the hlter. For the same reason, the 
noise-equivalent power of a phononic-isolated bolometer estimated in hgure |^b) is an 
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6. Conclusions 

We have explored the phonon hltering properties of aperiodic hnite-sized multi-moded 
phononic waveguides. An aperiodic hlter is composed of several cascading phononic 
stages of different size. The phononic geometry employed in each hlter stage results in 
complete bandgaps when simulated in the periodic limit. Thus, the hlter rejects acoustic 
and optical modes, and the synthesized phonon stopband may readily exceed 100% in 
fractional bandwidth. 

For the target application of bolometric sensors operating at 100 mK, an aperiodic 
phononic hlter with a stopband from 1.6-10.4 GHz reduces the thermal conductance of 
a uniform beam by a factor of 5 at 100 mK, and by an order of magnitude at 75 mK. To 
achieve a similar conductance with boundary-limited scattering requires millimeter-long 
and micron-wide beams, in contrast to the compact phononic hlter <10 /im in size. 

As envisaged, the aperiodic hlter geometry has self-similarity. While the fractal 
dimension is not a unique descriptor of the hlter prohle, it may be useful in future 
studies when relating the geometric complexity to the hltering performance. As a 
result of the multi-moded nature of the problem, merely stating the size ratio of the 
largest to the smallest phononic shape in an aperiodic hlter is not sufficient to describe 
its performance. Interestingly, if the geometry of the PnF changes according to a 
scaling law (e.g. logarithmically), the hlter transmission and held pattern radiated into 
a higher dimensional elastic space would remain constant over the desired frequency 
bandwidth 

A multi-stage coherent hlter can be applied in photonic [H], acoustic, and surface 
wave media. A wide stopband phonon hlter is particularly useful for thermoelectric 
generation, and for isolating high-Q mechanical resonators from thermal huctuations in 
cryogenic optomechanical devices. 
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Appendix 

We assume the following elastic constants for a parent material that roughly corresponds 
to silicon in the <100> direction [22]: E = 170 GPa, p = 2330 kg m“^, and z/ = 0.28. 
All unconstrained boundaries are set as stress free. The hnite-element computations 
inherently include coupling to all modes supported by a structure, including evanescent 
modes. 
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Since the displacement and traction vector fields are complex valued, care must be 
taken to extract the phonon transmission coefficient for a given solution in a finite- 
element computation, i.e. the location of an excitation source matters relative to 
scattering regions such as step changes in geometry. Edge loads are used to excite 
the four acoustic modes. To investigate the level of interaction between the edge source 
and reflected fields, a uniform beam with perfectly matched layers at both ends (denoted 
by AA for absorbing boundary conditions) is first simulated to “calibrate” the response 
of the beam to a given load. The simulation is repeated with one end of the beam fixed 
{AF). The reflection coefficient from the fixed boundary is evaluated from the complex 
displacement field, r = Uaf/Uaa — 1, and its deviation from unity is taken as a measure 
of the interaction strength between the source and the reflected fields. This procedure is 
generally known as a Thru-Reflect calibration [23] , and suggests an interaction strength 
of -30 dB at most between the source and reflected torsional acoustic mode power, and 
less than -50 dB for the other three acoustic modes. 
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